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Abstract 



We present a full study of the 3-body problem in gravity in flat (2-|-l)-dimensional space- 
time, and in the nonrelativistic limit of small velocities. We provide an explicit form of 
the ADM Hamiltonian in a regular coordinate system and we set up all the ingredients 
Oh| for canonical quantization. We emphasize the role of a U{2) symmetry under which the 

Hamiltonian is invariant and which should generalize to a U{N — 1) symmetry for N bodies. 
This symmetry seems to stem from a braid group structure in the operations of looping of 
particles around each other, and guarantees the single-valuedness of the Hamiltonian. Its 
role for the construction of single-valued energy eigenfunctions is also discussed. 
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^ ! 1 Introduction 
rn 

. The gravitational problem with or without matter in 2+1 dimensions [Ij has received consider- 

able attention in the past years (for reviews, see Ref. [2l|3]) as a laboratory for a nonperturbative 
treatment of gravity at the classical and quantum level. Recently, much work has been devoted 
to topologically massive 2-|-l-gravity in AdS spacetimes (see e.g. [1] and references therein), 
but several aspects of the more conventional problem with matter in an open spacetime are still 
^ ' unsolved, and addressed here. 

■ In fact, despite several efforts, an explicit quantum mechanical treatment of 0(2, 1) gravity 

with matter has been found only in the 2-body case, following the Deser-Jackiw-'t Hooft (DJH) 
classical solution [5j and its quantization [6l [7j . Subsequently, progress in the treatment of the 
classical A^-body case with a regular metric has been achieved both in the first-order formalism 
[SI [9] and, more implicitly, in the canonical one [H] by using a York-type gauge. But no 
conclusive work on the solutions of the A^-body quantum gravity problem is yet available for 
iV > 3. 

The purpose of the present paper is to perform further steps in the direction of a canonical 
quantum treatment of the = 3 case. By using the simplifying assumption of small velocities 
[101 [9] , we are able to provide the (single- valued) Hamiltonian and a complete set of constants 
of motion in a fully explicit form and in a regular coordinate system. Canonical quantization is 
then, in principle, straightforward. However, we have not been able, so far, to implement the 
monodromy condition on the energy eigenfunctions, and thus to construct the canonical Hilbert 
space. 

We use, throughout the paper, a regular coordinate system, such that the metric be single- 
valued everywhere, including in the neighbourhoods of (pointlike) external matter. On the other 
hand, it is known [1] that space-time is flat outside matter sources, so that one can use instead 
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Cartesian coordinates - characterized by various deficit angles - in which the conjugate momenta 
to the particles' positions are constants of motion. The mapping from regular coordinates in a 
York-type gauge to Cartesian coordinates was constructed in Ref. |9| thus providing an explicit 
expression of the constants of motion both for N = 2 and in the nonrelativistic A'^ = 3 case. 

The canonical formalism was then set up in Ref. [11] by providing the form of the Hamiltonian 
for N = 2 and its implicit definition for > 3. Its general interpretation is that of the scale 
factor in the asymptotic Liouville field occuring in the ADM parametrization of the metric. 
The external masses act as sources of the Liouville field, together with the so-called apparent 
singularities [12] of the problem. Therefore, the asymptotic scale factor is a function of the 
particle masses and coordinates which can in principle be computed. 

Our first task here is to provide the explicit form of the Hamiltonian for the 3-body case, on 
the basis of its definition in Ref. |llj . We are able to do that in the small- velocity limit, where we 
show that the Hamiltonian is simply related to a sum of squared moduli of two properly-defined 
relative momenta (called P3 and P2 in the paper), whose expressions in terms of regular particle 
coordinates and momenta are explicitly given in Sec. El Since such expressions carry branch 
cuts, there are nontrivial monodromy transformations for P3 and P2 when the particle positions 
(of particles ^2 and ^^3, say) turn around each other. Nevertheless, the Hamiltonian is left 
invariant (and is thus single-valued) because it possesses a U{2) symmetry and, moreover, the 
pair (^3,^2) turns out to transform as a U{2) spinor. 

The U (2) invariance of our problem has thus an important role in assessing the monodromy 
of the Hamiltonian, and the exchange symmetry of its equation of motion. We feel it should 
have a role in quantization as well, because it regulates the degeneracy of the wave functions 
and their monodromy properties. 

We are thus able to provide explicit solutions of the classical Hamilton equations and to 
express them in terms of Cartesian coordinates and momenta. Such solutions agree with those 
found in Ref. [9] in the first-order formalism. Given such explicit understanding of the Hamilto- 
nian structure, canonical quantization is in principle straightforward, and is here formulated by 
using a proper ordering prescription for the expression of the Hamiltonian in terms of regular 
coordinates and momenta. However, while the Hamiltonian is single valued, its eigenfunctions 
- characterized by a large degeneracy - are generically multiple-valued under the braiding of 
the 2 and 3 labels, for instance. We have some ideas on how to possibly construct monodromic 
eigenstates, but we have not found a successful procedure yet. 

After summarizing previous work on the canonical formalism and describing how to calculate 
the Hamiltonian in the nonrelativistic limit in Sec. [21 we provide its explicit form for = 3 in 
Sec. El where we discuss its U{2) symmetry also. The equations of motion and their classical 
solutions are given in Sec. HI the comparison with previous calculations in a different formalism 
is done in Sec. El while in Sec. El we sketch the quantization of the problem. We summarize 
our results and outline a few ideas which may prove useful in order to construct the quantum 
Hilbert space in the final section. Some technical details are discussed in the appendix. 

2 Hamiltonian formulation and nonrelativistic limit 

We want to describe the motion of N pointlike massive particles in a 2+1 dimensional open 
universe in a Hamiltonian formalism, in order to be able to perform canonical quantization. 
The complete formalism being available in several places in the literature \13\ [TT] (see also 
[HUH]), we shall not enter the details of the derivation. We provide the essential formulae in the 
general case in the first subsection, and specialize to the nonrelativistic limit in the next one. 
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2.1 General Hamiltonian for the motion of pointlike massive particles 

Let us start with the action in the Arnowitt-Deser-Misner (ADM) formulation |15^ [T6] . Parametriz- 
ing the generic hne element ds in the standard form 



ds^ = -N^dr + hijidx' + N'dt){dx^ + Wdt), (1) 



the action reads 

S 



where qn is the position of particle number n (p„ is the conjugate momentum), 11 is the canonical 
momentum conjugate to the spatial metric /i, and 



n 



IZ is the intrinsic curvature of the 2-spacelike slice and V is the covariant derivative compatible 
with the 2-dimensional metric hij. k? is related to Newton's constant through = SvrG, and 
has dimensions of an inverse mass in 2-dimensional space. It was set to ^ in Ref. |131 [TT] , and 
to 1 in Ref. (Sj [U [9] . We refer the reader to [171 ttS] for a complete expression including the 
boundary terms Sb- 

We can check that the variation of the action in Eq. ([2]) with respect to the metric components 
N , N^, hij and to the momentum Hij yields the Einstein equations in a first-order form. In 
particular, one gets the equations 71 = and Hi = which are the so-called Hamiltonian and 
momentum constraints respectively, and which we shall analyze below. The Hamiltonian for the 
particles' motion is obtained after all variables have been expressed as a function of and qn. 

We choose the York instantaneous gauge in which the intrinsic curvature K is zero every- 
where, which has proved useful in this context. Introducing complex notation z = + ix"^, 
Zn ^ qn + iqn: ^ud Pn = {p\ — ip^)/2, we impose the further gauge-fixing conditions: h^z = 
and hzz = 0. The spatial line element then takes the conformally flat form 

dt^ = e^"\dz\'^ (or, equivalently, hij = e^"5ij). (4) 

This equation defines the field a. 

The momentum and Hamiltonian constraints can now be solved for H and a. The momentum 
constraint Tii = determines the components of the energy-momentum tensor: 

27r ^ z- Zn 

(up to an entire function of the z variable which has to be set to zero to ensure good asymptotic 
properties of the metric, see Ref. [I3]). We will work in a center-of-mass frame in which YlnP"- ~ 
0. In this case, H is the ratio of two polynomials of respective degree N — 2 and N. Let us 
denote by za the zeros of the polynomial in the numerator. They are functions of the canonical 
variables Zn and p„. We write 

2vrn„(^-^n)' 
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where the normaUzation D is the so-called dilation factor 



-D = ^ ZnPn, (7) 
n 

whose imaginary part is half of the total angular momentum of the system of particles. 
The Hamiltonian constraint Ti = reduces to the nonlinear equation 

A{2a) = -2|2K2n|2e-2<^-47r^^„(52(z-z„), (8) 

n 

where A = d'^ + dy = ^dzdz and where we have introduced the dimensionless masses = 
K^m„/27r. Defining 

e"^* = 2|2K;2n|2e-2<^, (9) 
the new field a obeys a Liouville equation: 

A(2a) = -e-2^-47r^(/i„-l)(52(z-z„)-47r^52(^_^^)_ ^^^^ 

n A 

There are N sources located at the particles' positions Zn- The remaining sources za are so-called 
apparent singularities: They stem from the zeros of the momentum tensor 11. Their number is 
A'^ — 2 in the center-of-mass frame. 

The Euler characteristic of the 2-surface is an important parameter. It reads 

27r^ = - / (fzVhn = -- I (fzA{2(T) = lim -- /v(2cj) • dn. (11) 

2 J 2 J r^oo 2 J J. 

The notation for the integration element is d'^z = dxdy = dzdz/{2i). The last equality, obtained 
from the Gauss law on a circle of radius r with exterior normal n, provides the asymptotic 
behavior of a in the form 

e ~ -- , (12) 

kl»|A| A 

where A is a distance scale. 

The parameter /x is a constant of motion and is interpreted as the (rescaled) total mass of 
the universe. If all particles are static, there is no interaction energy, since in 2-1-1 dimensional 
gravity, spacetime is flat outside the sources, and thus there are no local interactions. Hence in 
the static case, the total mass of the universe is the sum of the masses of all particles: /x = Ylm 

It was found in Ref. [llj that the physical Hamiltonian is related to the logarithm of the 
scale A of the spatial distances, namely 



^=^ln|Ap (13) 

up to a constant term, where A is a function of the canonical variables which are the particle 
momenta p„ and positions Zn- To give sense to Eq. (jl3p . one has to make A dimensionless by 
dividing out a basic length scale (which does not appear to be fundamental, but rather related 
to the initial conditions.) In the following, all length variables (A, z, Zn- ■ ■ ) will be considered 
dimensionless. 

The Hamiltonian for the 2-body problem can be computed exactly from Eq. (jl3p . How- 
ever, for larger values of A'", the exact resolution of the Hamiltonian constraint ^ is already a 
formidable task |18] . Therefore, in this paper, we shall stick to the nonrelativistic approximation 
which was proposed in Ref. [9], and which we shall now derive in the Hamiltonian formulation. 
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2.2 Nonrelativistic limit 



To characterize the quasi-static approximation, we have at our disposal the + 1 parameters 
/i and fin- When the particles are at rest, the difference between the total mass of the universe 
and the sum of the masses of the particles 



f^-^f^n (14) 



is zero. So in the following, we shall always restrict ourselves to the lowest order in e. 

Let us start from the definition of the Euler characteristic (llip . Inserting Eq. ([8]) in Eq. (llip . 
we obtain the following relation for a defined in Eq. ([9]): 

— f <fze-^^ =e. (15) 
2vr J 

An analysis of Eq. (jlOp shows that on the particle singularities, the term e"^*^ has to behave 
like \z — Zn\'^'^^"~^^ ■ (This is consistent with the fact that e^^*^ vanish for z = Zn, as long as the 
masses /i„ are positive). In the right-hand side of Eq. (jlOp . the (5-functions then dominate at 
the particle singularities. Choosing a regularization, we see that from Eq. (llSp . the term e"^"^ 
is of order e. The scale e disappears when one takes the Laplacian of its logarithm as in the 
left-hand side. So setting e to zero, we see that the nonlinear term drops out and the Liouville 
equation for a boils down to a Poisson equation. The solution of the latter reads 

N N-2 
71=1 A=l 

where K \s a, complex number independent of z. We note that this solution is consistent with 
the one found in Ref. There, the analysis was performed in a first-order formalism, using the 
scale of the Cartesian particle velocities as a small parameter. We shall make the comparison 
sharper later on. 

To obtain the Hamiltonian, we only need the large-distance behavior of e~'^^ which can be 
deduced from the behavior of e^^*^ and of 11 using Eq. ([9]). First, from Eq. (I16p . 



e - ~ (17) 



-25- 

r 

\z\-^oo 



Second, from Eq. ([6]), the momentum 11 behaves like D/z"^, where D is defined in Eq. ([7]). 
Identifying the asymptotic behavior of e~'^^ just found with Eq. (|12p at the lowest order in e, 
we get an expression for |Ap from which, with the help of Eq. (|13p . we deduce the following 
formula for the Hamiltonian: 

up to an irrelevant constant. The parameter \K\^ can be determined as a function of the 
canonical variables and of the constant e using Eq. (jl5p . So our task is now to integrate Eq. (jl6p 
over the whole complex plane. An appropriate change of variable enables us to cast the integral 
in the form 

^i-li^|2|z2i|2M-2 j - ep^^-'ica - i?^'-^ ■■■\Cn- ci''^--' n - ^1' (19) 



A 



where we have defined the new variable Qn/A = {^n/A~'^i) / {^2—^1)^ -^ij = Zi—Zj and /i = Ylin=i 
consistently with the quasi-static approximation. 
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It is easy to check that these formulae ahow to recover, up to an additive constant, the 
Hamiltonian of the 2-body problem, which was written down by several groups [5t[6|[7t[8l[9|lll]: 

H = ^{\n\2K^p\^ + ^ji\ii\z\^) , (20) 

where p is the momentum p2 of particle #2 (keeping in mind that P1+P2 = 0) , z = Z21 = Z2 — Z1. 



3 Explicit Hamiltonian for the three-body problem 

So far, we have obtained a general formula for the Hamiltonian in the nonrelativistic limit 
as a function of the canonical positions and momenta of the particles, described in a regular 
coordinate system by the variables z^ and pn respectively. It is given in Eq. (|18|) . with D defined 
in Eq. ([7]) and \K\'^ obtained from the evaluation of the integral in (jl9p . In the present and 
the following sections, we shall specialize to = 3 bodies. In this case, the nonrelativistic 
Hamiltonian may still be expressed in terms of known functions and, therefore, may be studied 
completely. 

Interesting new features appear with respect to the two-body case. In particular, the Hamil- 
tonian possesses a U (2) invariance, related to a braid group structure of the particle exchanges 
and loopings. 



3.1 Explicit calculation 

Let us first gather the ingredients for the computation of the Hamiltonian (jl8l) . namely the 
expressions of \K\'^ and of D. 

We specialize Eq. (fT9]) to the case of 3 bodies. The rescaled position variables of particle #3 
and of the apparent singularity respectively read 

^ _ ^ ^31 > {P2+Pz)C, , . 

C = C3 = — , Ca = ■ -r- (21) 

We choose the origin of the frame at the position of particle #1 in such a way that zi = 0. Then 
Eq. ([19D reads 



(22) 



where /i = ^1 + /i2 + in the right-hand side of this equation. The dilation operator defined in 
Eq. dZD is 

D = Z2{P2+P3C)- (23) 

We now compute the integral that defines e (Eq. (j22p ). To this aim, we first expand 

le - CaP = - CP + (C - C)(C - Ca) + (e - C)(C - Ca) + |C - CaP (24) 
in order to be able to cast e in the form of a sum of standard integrals: 

e = |Kp|z2p^~' [Jii + (Ca - 0^10 + (Ca - O^oi + ICa - Cl'^00] , (25) 

where 

J^^ = y"d2^|^|2(Mi-l)|l _^|2(M2-1)(-^_^)M3-1+5(-^_^)M3-1+5_ (26) 
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This kind of integrals appears in conformal field theory and were computed e.g. in Refs. [201121]. 
The J- functions can be expressed with the help of hypergeometric functions: 



Sl3 

+ (-l)'^"'^5^,2,^^+^3+5_iF(l - ^3 - 5, 2 - ^1 - /i2 - /is - 5, 2 - /ii - /is - 5, C) 

^ ^/.2,Mi+M3+5-l^(^ - /i3 - ^, 2 - /Xi - /i2 - - 5, 2 - /ii - ^^3 - ^, C)^^, (27) 

where Si = sinvr/ij, Sij... = sin7r(/Zj + + ■ ■ ■) and -B«,/3 is the standard Euler beta function. 
Note that since we are working at the lowest order in e, we may always replace ^ /x, by the 
total mass ^ in each term that appears in the r.h.s. of Eq. ()25p . We define the four following 
functions: 

ha = B^^^^j,.^C,^^^>^''i~^F{^xi,l - ^2,^J'l+ hb = ha\^JLz^l+^JLz, ^23) 

/2a = -B^i+/i3-l,At2-^(l - /X3, 2 - /i, 2 - /ii - /i3, C), f2b = /2aU3^1+^t3- 

The functions /3a and /2a of the variable Q are two independent solutions of the hypergeometric 
equation 

C(l - Of'L + [2 - /^l - /i3 - (3 - /i3 - Mia - (2 - - ^l^)fia = (29) 

and, similarly, f^^ and /2fe solve the hypergeometric equation obtained from the previous one 
after having performed the shift ^^3 — )• ^^3 + 1. 

After some algebra, involving in particular the trigonometric identity 

sin vra sin 7r(Q + /3 + 7) = sin 7r(a + /3) sin 7r(a + 7) — sin 7r/3 sin 7r7 (30) 

with a = Hi, /3 = fi2 and 7 = /i3, we arrive at an expression of e/\K\'^ in terms of the /'s: 
4g 

where 



k2p'^-' [iVsKCA - C)/3a + /36P + A^2|(Ca " C)/2a " f2b\^] , (31) 



The Hamiltonian (1181) then reads 



iV2 ^ iV3 ^ (32) 

7rSi23 TTSlS 

^ = ^1 In + (/X - 1) In k2|' + In [Ns\{Ca - Ofsa + /sfel' + A^2|(a " C)/2a - f2b\^] }. 

(33) 

Introducing some more notations, whose physical interpretation will be given below, the Hamil- 
tonian (j33p may be rewritten in a compact form 



1 (2.CY(|P2P + |P3P) 



where we have defined 



P2 = Dz^'WN2 [f2b - (Ca - C)/2a] and P3 = Dzi^-'^jN^ [h, + (Ca - C)/3a] (35) 
or, written in terms of the canonical variables Z2,C and P2,P3 with the help of Eqs. (j2ip and (|23p . 

-P2(2:2,C,P2,P3) = z!^F2{C,P2,P3), F2 = a(C)P2 + /3(C)P3 ^^^^ 
-P3(^2, C,P2,P3) = ^2-^3(C,P2,P3), ^^3 = l{C,)P2 + SiQps, 
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the coefficients a, 13,^,5 in front of the momenta being defined as 

a = \/iV2/26, /3 = V^C(/26 - (1 - Of 2a), 



(37) 

7 = ^/Nzhb: S = VN3C{f3b + (1 - C)/3a). 

This set of notations will help to write the calculation of Sec. |3]in a simpler way. 
3.2 Interpretation of P3 and P2 in the small mass limit 



Before embarking with the study of the Hamiltonian (j34p . we wish to try and give an interpre- 
tation of its expression. It is quite straightforward in the small mass limit ^„ ^ 1 for n = 1, 2, 3 
in which gravity effects vanish. Let us write the expression of P2 and P3 in Eq. ()36p in this limit. 
We start with the /'s defined in Eq. (|28p : At the lowest order in all the they boil down to 

ha — T, J2b — , /3a — 7 H r, J3b — , (oSj 

and it follows that 

p _ 1 ^2/^13 _ Pl±P3\ p 1 ^^1^3 (P^ _ P}\ .2g^ 

^ ^ \/^ /i Atl3 / ^ ^ \/^ /^13 VA^3 /^l/ ' 



Thus we see that P3 is the relative momentum of particles #1 and #3 with the normaliza- 
tion 1//V3 ~ \/ /Ui/i3///i3, and P2 is the relative momentum of particle #2 with respect to the 
(13) subsystem, with the normalization \//V2 — ■\/a^2A*13/V- the absence of interaction, the 
nonrelativistic kinetic energy of 3 pointlike particles reads 

p = 2K^f|L + A + J_U4.^fM + N!^N!y (40) 

V2/ii 2/X2 2/^3/ \ /ii /12 /i3 / 

With the help of the expressions for P2 and P3 in the small mass limit given by Eq. (|39p . it can 
be rewritten as 

P = 4k2(|P2|2 + |P3|2). (41) 
Note that this kinetic energy does not coincide with the Hamiltonian in Eq. (|34p . which reads 

H = ±-\n^. (42) 
2k2 4e ^ ' 



This is actually related to the time gauge which was chosen, and is better discussed in Sec. 14. 3[ 
3.3 Monodromy properties, U{2) symmetry and braids 

We want to check that the Hamiltonian (j33p is well-defined. It is necessary that it be single- 
valued, i.e. invariant under the loopings of C, around the branch points at 0, 1 and cxd, and 
invariant under the exchange of the labels of the particles. 

3.3.1 Monodromies and U{2) symmetry 
Let us introduce the two objects 

in such a way that the momenta (P3, P2) defined in Eq. ([35]) can be conveniently rewritten as 

P = {^^^ = Dz^-' K + (Ca - CW] ■ (44) 
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We compute the monodromy matrices M^2 and M31 which correspond to the loopings of parti- 
cle 7^3 around particle ^2, and #3 around 7^1 respectively. These transformations amount to 
substituting Q — 1 ^ e^*'^(C — 1) and C, — )■ e^*'^^. We find that Ua and ab transform according to 
the matrix (see Appendix Rl for the details) 

for the looping of particle #3 around particle ^^1. As for the looping of particle #3 around 
particle #2, the corresponding transformation reads 

M32 = e*-(^2+'^«) f ^32\ ^.^j^ fa32=COS7r(/i2 + M3)+^sin7r(A.2 + /^3) cosa32, ^^g^ 

V-^32 032/ |632 = zsin7r(/i2 + ^3) sina32, 

and where the angle 032 has been defined as 

S1S2 - S3S123 . 

cosa32 = , sma32 = . (47) 

S13S23 Sl3'S23 

The matrices M31 = M3ie"*''(''i+^3) and M32 = M32e-*''(^2+M3) are SU{2) transformations, so 
that the monodromy group generated by M31 and M32 is actually a subgroup of {7(2). 

The object [Qa — C)(^a also transforms according to M31 and M32, and the same holds for 
the sum [Qa — C)'^a + ^b- This implies that the norm of this SU{2) spinor, which reads 

A^3|(Ca - Oha + hb? + N2\{U - Ofla " hb? (48) 



is stable under the monodromy transformations. Consequently, the Hamiltonian (|33p is single 
valued under all possible loopings of particle #3 around the two other particles. 

We have not considered the looping of particle #2 around particle #1. The corresponding 
transformation of the spinors can be computed: The best is to exchange #3 and #2 in such a 
way that the looping of particle #2 around #1 be the substitution Q — s- e^*'^^, and eventually, to 
go back to the initial coordinates by applying the inverse exchange transformation. Therefore, 
we postpone this discussion after we have studied how the momenta change when the particles 
are relabeled. 

3.3.2 Relabeling symmetry 

Exchanging the labels of the particles yields nontrivial transformations of the spinor (i-*3,P2)- 
We now check that the Hamiltonian is invariant under these transformations. 

We first perform the exchange o //3 and {zi^pi) -H- (-23, P3) in the expression of the P's 
in Eq. (|36p . Since the details are lengthy, we defer the full calculation to Appendix |Al We 
find that the transformation of the spinor (^3,^2) can be written as the multiplication by the 
diagonal matrix 

^31 - e Q g-i7r(Mi+M3)/2 j ■ ^'^^l 

In the same way, the exchange of the labels of particles #3 and #2 is represented by the 
multiplication by the matrix 

MA«2(m) = e^"(^^+^«)/2 /g-i^(Ml+M2/2+M3/2)^^ e^'^(^'3-M2)/2^^^ 

"^^^ " ^/^T^ V e^'^(^2-M3)/2^/5YiI^ _e^-(Mi+/^2/2+A.3/2)^/^ 



We have written the two matrices as the product of a U{1) phase and a U{2) matrix of deter- 
minant — 1. (This is clear in the former case, and can be easily checked using the trigonometric 
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identity (j30p in the latter case.) It is then obvious that the spinorial norm + is left 
invariant by these transformations, and since e is trivially invariant, the same is true for the full 
Hamiltonian ()34p . 

We check that the monodromy transformations are the "squared" of the relabeling transfor- 
mations. More precisely, 

^sT'^sT' = ^31 and r3^|^3(Mi)^A«|M2(Mi) ^ (5^) 

We are now in a position to easily get the last monodromy transformation which we left 
uncomputed in the last subsection, namely the action of the looping of particle #2 around 
particle #1 on the spinor (P3, P2) that we shall denote by M21. We write 

IVI21 — [T^2 ) ^31 ^31 ^32 • W^J 

The result of the matrix multiplication may be written in the same form as Eq. (j46p . except for 
the values of the U{1) phase and of the SU{2) angles: 

^ gMm+M2) ^.^j^ faai =cos7r(/ii+^2) + isin^(^i+^2) cosa2i, 

where 



S2S3-S1S123 . 2^SiS2S3Si23 f^,. 

cos 021 = , sma2i = ■ (.54) 

S12S13 S12S13 

The monodromy M21 is nothing but the exchange of particles ^2 and ^^1 repeated twice, which 
then has the following matrix representation: 



21 



/'^-1\M3M2(Mi) ^M2/tl^M3M2(Ml) /gg-v 



and whose explicit expression reads 



3.3.3 Braids 

There are interesting relations between the transformations (|45p . (|46 p , (|53 p . (|49 p . (|5Up . (|56 p which 
seem to point to some underlying braid group structure [22j. (It was anticipated by 't Hooft 
in Ref. [B] that braids and/or knots would play a role in this problem). We mention here the 
correspondence between the monodromy group and the braid group since it could be important 
for quantization, but without entering the details: We shall refer the interested reader to the 
litterature for an introduction to braids. 

Let us consider 3 strings, each of them being attached to one of the particles, and which 
carry the mass of the latter. We assign the relabeling matrices r to the crossings of two out of 
the three strands. For each crossing, there are two possibilities, depending on which one of the 
particles passes above. We choose the convention that the matrix the strand 

from the middle of the braid associated to particle #i of mass /ij above the strand initially at 
the bottom associated to particle ^^j. The matrix tI^^^'' takes the strand at the top above the 
strand initially in the middle. The opposite orderings are represented by the inverses of the r 
matrices, namely 

= (^M,M.(/.fe)^t and {t^^^T'^" = {ri^i^''')^. (57) 
Note that one has to keep proper track of the particle masses along the strands. 
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Now we check by an explicit calculation that the following relation hold: 



'31 



(58) 



This identity is reminiscent of the defining property of the -B3 group of braids on three strands, 
see Ref. [22]. 

Given the relation between monodromy and relabeling transformations, it is not difficult to 
check that the group generated multiplicatively by the matrices M31, M32 and M21 is homo- 
morphic to the pure braid group on three strings. The latter is the subgroup of the braid group 

which preserves the ordering of the strands. We check by explicit matrix multiplication that 
the two following defining relations of the pure braid group (see the corresponding equation in 
Ref. [22j, which differs only by an appropriate relabeling) are identically verified: 

M31M32M3-I = M2-/M32M21 , M31M21M3-/ = M2-/M32^M2iM32M2i. (59) 

4 Equations of motion and conservation laws 

So far, we have derived and studied the Hamiltonian which describes the evolution of three 
pointlike particles (see Eq. We are now going to investigate deeper the dynamics of the 

system. Prom the Hamilton-Jacobi equations, we will be able to compute the time evolution of 
the dilation factor (i.e. also of the total angular momentum) (Sec. 14. ip . We will show that the 
momenta P3 and P2 in terms of which the Hamiltonian is written are conserved (Sec. 14. 2p . and 
we will define the positions Z3 and Z2 canonically conjugate to P3 and P2 respectively (Sec. 14.3]) . 

In order to avoid to have to carry along 2/t^ factors, let us take the momenta dimensionless 
by setting 2 k? = 1. 

4.1 Hamilton-Jacobi equations and time evolution of the dilation factor 

In our definition of complex positions and momenta (see Sec. [2]), the Hamilton-Jacobi equations 
for the time evolution of the coordinates read 

dH dH 
Z2 = ^— , = — , (60) 
dp2 dp3 



that is, from Eqs. ([MD and (l36]l : 

^ aP2 + 7^3 . ^ m + 6P3 

^2 = ^2 |p^|2^|P3|2 ^ ^3 = ^2 |p^|2^|P3|2 - (^1) 



Comparing to the definitions (j36j) . one immediately sees that 

P2Z2 + P3Z3 = 'i^- (62) 
Similarly, the Hamilton-Jacobi equations for the evolution of the momenta read 

dH . dH 

P2 = P3 = (63) 

OZ2 OZ3 



Replacing H by its expression ()34p . one gets, after the further replacements of P2 and P3 by 
Eq. dMI): 

■ ^2 fdF2^ , dF3- 

''-~2- IP2P + IP3P te^^ + d^"^', ... 

- "IP2P + IP3P V 9^^^ + 9^^^ 
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Since F2 and F3 only depend on (, their derivatives with respect to Z2 and Z3 read 

dz2 Z2 dC ' dzs Z2 dC ' 
We then see that the following conservation law is satisfied: 

= -Ai- (66) 

Combining Eq. (|62p and (|66|) . we find that the dilation factor D defined in Eq. (|23p has a linear 
evolution with the time t: 

D = l-H. (67) 

4.2 Cartesian momenta 

We are going to prove that P2 and P3 are constants of motion. To this aim, we need to compute 
the Poisson brackets of the Hamiltonian with the momenta: 



4=2,3 

and similarly P3 = — P3}. We easily find 



p ^ {P2,P3}P3 p ^ {P3,P2}P2 . . 

I^2P + |P3P' ' |^'2P + |P3P' ^ ^ 

We are going to show that the Poisson bracket {P2,P3} vanishes. Replacing P2 and P3 by their 
expressions ([36]) . we get 



{P2, P3} = zl^^^ I {cl{8 - 7C) - 7(/3 - «C)] V2 

+ [/3'(<5 - 7C) -nfi- «C) - /i(«<5 - /37)] P3 \ . (70) 



From the definitions of a, /3, 7, 5 in Eq. (j37p . the following identities hold: 

/3-«C = -\/^C(l-C)/2a 

5-7C= V^C(l-C)/3a (71) 



a5-/37= Vi^C(l-C)(/2a/3fe + /2b/3a). 

The derivatives of a, /3, 7, 5 can be expressed as a function of the /'s and the second derivatives 
of hh and by using 

/26 = -[i-iha, fzh = fJ'3f3a (72) 

which are consequences of standard identities between hyper geometric functions. 
We insert Eqs. ([71]) , ([72]) into Eq. and get 



/^3 I 

+ /3a[C(l-C)/2; + /^3(l-/^)/2j }. (73) 

Thanks to the hyper geometric equation (p9]) applied to /3b and /2b, we see that the term under 
the square brackets cancels identically, and thus {P2,P3} = 0. Equations (f69]) eventually show 
that P2 and P3 are constants of motion: 

P2 = 0, P3 = 0. (74) 
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4.3 Cartesian positions 

We have just exhibited two constants of motion: The two relative Cartesian momenta of the 
particles. The Cartesian velocities, which are the time derivatives of the Cartesian positions Z2 
and Z3, should also be constant. We are going to define them as the variables conjugate to the 
momenta P2 and P3. 

The dilation factor D = P2Z2 may be expressed with the help of P2 and P3 in the form 

We define the new variables Z2 and Z3 in such a way that 

D = {l-fi){Z2P2 + ZsP3). (76) 

Let us work out the explicit expression for Z2 and Z3 as a function of the canonical variables. 
We introduce the notation 

W = ^/N^zihahb + /2 Jsa), (77) 

which is the inner product of the spinors a a and o"b, namely 

W = (c7b,a„) = a,^(^J 'q)"- C^S) 

and thus a SU{2) invariant. It turns out that W has a simple expression, see Eq. (5.5) in Ref. [9]: 
As a matter of fact, it is a Wronskian function for the solutions of a second-order differential 
equation. 

After comparison of Eqs. (j75p and (I76p . the explicit expressions for Z2 and Z3 are 



<'^^"^ (79) 

^'-'^ (i-/.)W 

The pair Z = {Z2, —Z3) makes up a spinor, in such a way that the spinorial product {P, Z) = 
P2Z2 + P3Z3 be an invariant under SU{2) transformations. 

In order to make contact with our intuition of classical physics, let us again take the small 
mass limit in which gravitational effects vanish. Then a straightforward calculation leads to 

Z3 = Vn'sZs, Z2 = (z2 - ^zs] , (80) 

that is to say, Z3 is the position of particle #3 up to a normalization, and Z2 the one of particle 
#2 with respect to the system of particles (13). (We recall that we have chosen a frame in which 
zi = 0). 

The dilation factor D has a constant time derivative and the P„ are constants of motion, 
we expect that the Z„ also have constant time derivatives. This turns out to be true, and 
furthermore, the Z^ s are canonically conjugate to PnS. The first point can be shown by 
evaluating 

7 [U 7 \ {^n,P2}P2 + {Zn,P3}P3 

Zn = -{B,Zn\ = |p^|2 _^ |p^|2 ' V«l) 

Let us perform the calculation completely for n = 2. The Poisson bracket of Z2 and P2 that 
appears in Eq. (jSTj) reads 

{Z2,P2} = zt^ [(1 - ix)aZ2 + (/? - aQZ'^ . (82) 
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dC \w J C(i-C) w 



To compute the derivative of Z, we use the more general formula 

d ( f'\ _ ^^3(1-/^) / ^ggsj 

valid for any solution / of the hyper geometric equation (j29p . We introduce the derivatives of 
/sfe in Eq. (i79]l expressed in Eq. (|72]) . One then uses the previous formula with / set to f^b to 
compute Z2 in Eq. ([82]) . With the help of the identities (j7T]) , one arrives at 

{Z2,P2} = 1. (84) 

Similar calculations lead to the following Poisson brackets: 

{Z2,Pz] = 0, {^3,^2} = 0, {Z3,P3} = 1, {^2,^3} = 0, {^2,^3} = 0. (85) 

Thus, 

= = \P2?^+\p.r ^^^^ 

For consistency, we easily check that D = \ — ^. 

It is possible to recover the usual Cartesian equations of motion 

by changing the time gauge as follows: 

^ = 2(|P2P + |P3P). (88) 

This reparametrization depends on the variables z„ and pn- With this choice for time and taking 
Zn and Pn as canonical variables, the T-evolution would be given by the Hamiltonian 

£; = 2(|P2|' + |^'3|'), (89) 

up to a constant. Such a Hamiltonian would be simpler and more intuitive since it is the 
nonrelativistic kinetic energy of the three particles in the center-of-mass frame (see Eq. (|4ip and 
the discussion above it), but the phase-space variables would not be single valued. 

5 Consistency with previous calculations and tentative exten- 
sion to many bodies 

In the present work, the computation of the Hamiltonian is based on a calculation of the total 
mass ^ = /i„ + £ as a function of the masses, positions and momenta of the particles in the 
framework of the second-order formalism. 

On the other hand, since (2+1) gravity is a topological theory, we know that pointlike parti- 
cles move on straight lines with constant velocities Vn, when appropriate (Cartesian) coordinates 
are chosen. The total mass ^ can be obtained by writing the total Cartesian momentum of the 
system of the 3 particles. The result for ^ should be the same as the one obtained in the previous 
section: This is what we are going to check here (Sec. 15. ip . This computation will also help us 
to establish the relationship between Vn and the derivatives of It can be extended to four 
(or more) particles and allow us to guess the form of the Hamiltonian in these cases (Sec. 15. 2p . 
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5.1 Comparison with the first-order formalism 

Following Ref. [9], we introduce a spin-^ representation of the holonomy related to transport on 
a curve around the particle n of mass /i^ and velocity Vn in the form 

r / T/ \ f^n bn\ , J a„ = cosvr^n + i7„sin7r//„, 
Ln{Hn,Vn)=[7 „ , where < • f> ■ v90) 

[bn = -z7„\/„sm7r^„, 

and 7n = 1 / y/ i-\Vn\'^ ■ In order to "measure" the total mass of the system, we may travel on a 
loop that goes around all the three particles. The total mass of the system ^ is then computed 
from the trace of the holonomy given by the product of the three nicitrices Ln- 

COSTTfl = ^Tv[L3{fi3,V3)L2{fi2,V2)Li{fli,Vi)]. (91) 

The order of the product is determined by the choice of ordering the particles anticlockwise in 
space. The calculation leads to 

cos TT// = cos TTfli COS TT fj,2 COS TT fl^ 

— (^7172(1 — Vi • ^2) sinvr/ii sin7r^2 cosvr/ia + [cyclic permutations] 

1 _ _ _ 

+ -7i7273(-^23Vi + V13V2 - V12V3) sinvr/ii sin7r^2 sin7r/X3. (92) 

We now expand to the lowest order in the velocities. We observe that the result can be written 
as the sum of the masses of the particles and of a quadratic form of the velocities, representing 
the total nonrelativistic kinetic energy of the system: 

/i = /il + /"2 + /X3 + l{Vl2 Vi3)Q (^^^^ , (93) 

where 



]_ ( S2S13 -e-*"'^^S2S3' 
2S123 V-e"''''S2S3 S12S3 

We may perform a standard decomposition of Q in a product of lower diagonal L, diagonal D 
and upper diagonal U matrices Q = LDU, with in this case, 

^=2(T ^J'^=(o ^ ^^^J andL = [/t (95) 

since Q is Hermitian. From this decomposition, we immediately write the kinetic energy as a 
sum of squared moduli: 



Z^^^n = -\ ^12 Vi3 + \Vi3\ 



(96) 



Note that the procedure used to arrive at this factorized expression may be repeated for a 
number of particles larger than 3. We sketch it in the next subsection 15.21 

A result of the first-order formalism developed in Ref. [9] is that the Cartesian velocities Vn 
can then be expressed with the help of the canonical positions z„ of the particles (see Eq. (4.19) 
in Ref. [9]). They read 

Vin = 2Kozl^-' / d^e^-H^ - IT'-HC - Cr'-\^ - Va), (97) 
Jo 
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where r]A is the apparent singularity. From the integral representation of the hypergeometric 
function 

dzz'^-^l- zf-^{l-tz)^-'^ = Ba^l3F{a,l--f,a + I3,t), (98) 



and after appropriate changes of variables, we express the V^s in terms of the /'s defined in 
Eq. ([28]), namely 

^12 = -2Koz^-'e-^-^- \ ^e-(^^+'^3) [^3^ + (^^ _ ^) J3J + [f^^ _ (^^ _ ^^f^^^ 

{ si3 J (99) 

Fi3 = -2Koz^-'e'^^'^-^-^^'^ [ht + (riA - Ofsa] ■ 
Inserting these expressions for the velocities in Eq. (j96p . we find 

f.-^f.n = 2\Ko\'\z2\'^-' (^\{rjA - Oha + h,? + ^|(r?A " C)/2a " ^ft^ (100) 

„ V^rsis 7rsi23 / 

Thus £ in Eq. pip matches — ^ /^i given by Eq. (jlOOp from the first-order formalism provided 

r}A = QA and |i^o|' = (101) 



Now we can establish the relation between the velocities, the Cartesian momenta and coordi- 
nates. Identifying Eq. ()99p with the expressions of the momenta (I35p . we write 



V21 = fe--i + e--.^^ 



On the other hand, by definition of the Cartesian momenta P2 and P3 



(102) 



\D\^ \P2\^ + \P3\^ 



(103) 



The right-hand side is a numerical constant since e has to be identified with the nonrelativistic 
kinetic energy in the Cartesian time gauge. Therefore 

^^31 = e'^ 



2P s 2P ^ 



N2 Sl3 ViVg 

where we have introduced the angle (p to absorb all irrelevant phases. One may now replace P2 
and P3 by the derivative of Z2 and Z3 with respect to T (see Eq. (187p ). Taking the small mass 
limit, with the help Eq. (j80p . it is easy to see that V31 and V21 coincide in this limit with the 
derivatives of Z3 and Z2 respectively. 



5.2 Guessing the form of the Hamiltonian for many bodies 

As was suggested before, we may express the total kinetic energy e of a system of many bodies 
as a function of the Cartesian momenta of specific subsystems of particles. Then from a change 
of the time gauge, we may infer the form of the Hamiltonian. We address in some detail the 
case of the four-body system. We write a relation similar to Eq. (I9ip : 

cos7r/i = ^Tr[L4(/X4,V4)L3(/i3, V3)L2(/i2, V'2)-Li(/ii, Vi)], (105) 
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which we then expand at lowest order in the Cartesian velocities Vu 

1 /^i^ 

M = ^1 + /U2 + ^3 + ^4 + -{Vl2 Vi3 Vu)Q h43 I • (106) 

In this case, the matrix Q encoding the quadratic form of the velocities reads 

Q = I -€'^"'^8283 S3S124 -6-^^^"^^ S3S4 . (107) 




S2S4 -e^'^^i^ggg^ S4S123 



Performing as in the three-body case a Q = LDU decomposition we find the following formula 
for e = ^ - ^^^^n■ 



S2S134 



-(— 

2 y7rsi234 



S134 -5134 



_^S3£l4 
7rSi34 



+^\n,n. (108) 



The identification with the classical equation for the total kinetic energy of the system (in the 
Cartesian time gauge) 

e = 2(|P2p + |P3p + |i^4n (109) 

leads (up to phases) to an expression of the Cartesian momenta P4, P3, P2 of particles number 
7^4 relatively to #1, ^3 relatively to the system (14), and 7^2 relatively to the system (134) 
respectively. This statement can be checked in the small mass limit in which the phases become 
irrelevant. We find for P4, P3 and P2 formulas similar to (j39p and which, up to normalization 
factors, are the expressions of the momenta of free systems of particles as a function of their 
Cartesian velocities. 

The strong similarity with the 3-body problem suggests that the Hamiltonian for = 4 
particles and more generally, for an arbitrary number of particles, would read 

H = in(!^!!±l^i!ll:::±lM), (no) 

Note however that the expressions of the P„ (or Vi„) in the regular phase space coordinates 
(zmPn) would involve line integrals in the complex plane with N cuts, which are not expected 
to be related in a simple way to known functions, as they are in the 2- and 3-body cases. 



6 Towards the quantum Hilbert space 

We impose canonical quantization by replacing the Poisson brackets {•, •} by the commutators 
—ih[-, •], which is realized by the substitution pj — )• —ihdzy In units h = c = 1, this results in 
the usual relations between the positions Zi and momenta pj, namely 

[zi,pj] = i6ij, [zi, Zj] = 0, [pi,Pj] = 0. (Ill) 

It turns out that the same rules apply to Zi and Pj , provided the ordering of Zi and pj is the one 
given in the defining equation (I36p (i.e. the position operators are to the left of the momentum 
operators). 
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The problem then arises of finding an exphcit expression for the eigenfunctions of the Hamil- 
tonian, so as to construct the quantum Hilbert space and the scattering amplitudes. It turns 
out that it is not difficult to find explicit solutions for the Schrodinger equation, because the 
Hamiltonian is a simple function of the P2, -P3 (and -P2, P3) operators, which are diagonalized 
by the "plane waves" 

where A;2, ^3 label the eigenvalues of -P3 and Z2, are expressed in terms of the regular 
coordinates Zj by Eq. (j79p . The corresponding energy eig envalue of Eq. (gSD is ^ = 2(|A:2p + 

However, the wave functions (jll2p are not single valued, but transform in a way induced 
by the monodromy transformations M of Sec. 13.31 In fact, if we let Z2 — >• = Z2e^^'^ or 
-23 — ^ = ^36^*'^ or combinations thereof, then (Z2, —Z^,) transforms as a spinor. The invariant 
combination in the exponent of (jll2p translates this transformation into a U{2) transformation 
of the momenta: k = {k^, k2) — )• k^'^ = {k^^ , ^2^)- The transformed wave function reads 

This wave function possesses the same energy eigenvalue E = 2(|fc2p + 1^3^) as '>pk{zj) in 
Eq. (I112p . In other words, ipki^j) transforms into another wave function of the same degeneracy 
class according to a U{2) representation of the monodromy. 

A similar situation arises already at the level of the simple two-body case, where the analog 
of Eq. (I112p is the wave function 

= e^['=^W+c-c.] ^ gi[fc.i-M/(i_^)+c.c.] ^ ^ ^ 2|A:|2. (114) 
Similarly, it is not monodromic, but transforms by a f7(l) monodromy as 

ipk{z) ^pkie^'^'z) = V'fce2«-«(^) . (a = 1 - /^), (115) 

which remains in the same degeneracy class. One may think that in the two-body U{1) case, it 
would be possible to construct linear combinations of plane waves which are single valued and 
have the same energy by just summing with unit weight over the whole monodromy group, as 
follows: 

+00 

Mz)= V'fc(e'^™^). (116) 

n=~oo 

However, this turns out to be too naive. Unless the monodromy group is finite (for fractional 
values of a), the series sums up to zero, after appropriate regularization. 

A solution which is both finite and monodromic was proposed by Deser, Jackiw [7J and 't 
Hooft [6j. It also consists in writing a superposition of plane waves degenerate in energy (which 
differ by the azimuthal angle of their momenta), but with nontrivial weights. Let us introduce 
some useful notations: 

X = \k\\z\'^/a , (3 = arg{k) , e = -argz. (117) 

Then 

Mz) = [ ^/(/3)e-™«(/^-^) (118) 

solves the Schrodinger equation when the contour C^g is a combination of Schlafli contours in 
the upper/lower complex plane of the (3 variable, namely the union of the broken line 

C_|_ =] — vr -|- ioo, — TT, TT, vr -|- ioo[ (119) 
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translated along the real axis by a9 and of its complex conjugate C_ (we introduce the notation 
C = C+ + C- for the union), with appropriate deformations in order to avoid the possible 
singularities of /. The function / is meromorphic and well behaved for Im(/3) — t- ±00. The 
convergence lines in /3 are f3 — a9 ^ (2n + l)7r it ioo, thus the contour Cae may be translated by 
I3± 27m. 

For each choice of C one may now construct a monodromic solution by replacing 9 hy 9n = 
9 + 2mT and summing over n: 



cos(/3-ae„) ^ / ^ 

c 2vr 



^/(/3 + «^n) 



The function 



= ^/(/3 + 27ran) (121) 

n 

which appears under the square brackets is periodic of period 27ra, so that ^/^^ is monodromic 
under the rotation 9 ^ 9 + 2tt. 

Such a monodromic expression can then be projected on the integer angular momentum m 
as follows: 



-irnd 



gix' cos /3 



(122) 

The integral over 9 (under the square brackets) yields the Fourier coefficient fm of the periodic 
function / times the phase e*™'^/". The integral over /3 is then seen to yield the Bessel function 
J\m\ja{x)-, which is the DJH [TJ |6] result. The role of the weight factor / is to suppress the 
contributions of the secondary sheets so as to avoid too much destructive interference. For the 
DJH scattering process / ~ l/(/3 — vra) and / ~ tan (/3/2a). 

Is a similar procedure available for U (2) (or more generally U{N — 1))? On one hand, when 
the monodromy group is finite, we may write 

Uz,) = E V'.m(Z(z,)) = E e^t'^f ^^(^^)+'=3-^3(..)+c.c.]^ (123) 
monodromies M M 

where k is now a U{2) spinor with components (^3,^2), and Z has components {Z2, —Z^). But 
on the other hand, despite various attempts, we have been unable to extend to this non-abelian 
case the general harmonic analysis of Eq. (I122p . The main obstacle is the nonabelian group 
multiplicative structure which does not allow an explict evaluation (with identification of good 
quantum numbers). 

Let us however give an example that shows how a simple structure may arise for proper 
(rational) values of the masses. Let us consider the case in which the particle masses are 

= /j,2 = 1^3 = \- Then the three basic monodromy matrices are related to the Pauli matrices 
cJi, (72 and (T3 through 

M31 = -CT3, M32 = ai, M21 = (72. (124) 

The group generated by the monodromies is finite in this case. It is made of the sixteen ^7(2) 
matrices el, eai, e(J2, ecr3, where e G {l,i, —1, — i}. 

In this particular case, the Cartesian coordinates have a simple expression as a function of 
the regular coordinates Z2 and C = 23/2:2 • Equations (f79l) indeed boil down to 

47r 1/4/. r. r\V2 „ 47r 1/4/. r. r\l/2 
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We now apply Eq. (jl23p to obtain the explicit monodromic wave function: 



ee{l,i,-l,-i} 



+ e 
+ e 



i[.fc24/*(i+^W)'/'-.fc34^'(i-^W)'^'+c.c.] 



^ gi[-iefc24/^(l-VT^)'/'+iefc34^'(l + V^)'^'+C.C.] I ^ (126) 

7 Summary and suggestions 

In this paper, we have fully understood the canonical Hamiltonian structure of the 3-body 
problem in the nonrelativistic limit in which the particle velocities are small. Let us summarize 
our findings, repeating the explicit formulas to which we have arrived. 
We have provided the expression of the Hamiltonian in the form 

H = —ln^- ^ ^' ' " ^ 127 

2k2 4e ^ ' 

where P2 and P3 are properly-defined relative momenta, which are given as explicit functions of 
the canonical coordinates {z2,P2) and (^^3,^3), while e is a fixed parameter which does not play 
a role in the dynamics: It represents the difference between the total (adimensionalized) mass 
of the Universe fi, and the sum of the masses /x„ of the particles, which are all constants. The 
nonrelativistic limit implies that it is small compared to 

By putting Eqs. (j36|) and (f37|) together, we obtain the expressions 



P3 = 4'0V3 [/36(P2 + CPs) + /3aC(l " OPs] , 
P2 = [f2b{p2 + CPs) - /2aC(l - Ops] , 

where the normalization factors and are given by Eq. (I32p . namely 

sin vr/ii sin 7r/i3 sinvr^s sinvr^ig 

-'V3 = : , iV2 = : , U^yj 

7rsin7r^i3 vrsmTr^ 

where fj-ij = Hi + fij and fJ- = fJ-i + fJ-2 + P-s- The /'s are solutions of the hypergeometric equation 
with specific coefficients (see Eq. ([29]) ). and their expressions read ([28|) 

r(/ii)r(^3) 

/3a = — p7 — ^ — C^'' Fifj-iA- fJ'2,mX), 

^^^^^ (130) 

r 1 - fl T{fil3)T{fl2) 00 /-\ 

f2a = WT-. ^(1 - M3, 2 - /i, 2 - /ii3, C), 

1 - /ii3 r(^) 

and /3f,, f2b are obtained from f^a, Ha respectively by shifting /U3 to /U3 + 1. P3 and P2 are 
constants of motion which are related to the Cartesian momenta. 

One can also build the Cartesian coordinates, which have relatively simple expressions. From 
Eq. ([79]), we get 



^^^^ -F(^3,/x- 1,^3,0, 



1 - /i \/i^r(^i3)r(^2) 

^^^'^^ -C^-'^-F(m2,1-W,2-/.i3,C). 



(131) 



i-ms V^r(/xi)r(/i3 
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Finally, the Cartesian time T is related to the ADM time t by 



dt 2(2k2)2(|P2|2 + |P3|2)' ^-'"^^^ 

which represents a time-gauge change. 

An important role is played by the U{'2) symmetry of under which (P3, P2) and (Z2, —^3) 
transform as spinors. This symmetry regulates the monodromies and in particular implies that 
H is invariant, that is, it is single valued. It is also useful to represent the exchange symmetries 
of the problem, e.g. the 2 -f^ 3 exchange, under which the Hamiltonian is also invariant. We have 
argued that the U{2) symmetry will be replaced by a U{N — 1) symmetry in the nonrelativistic 
A^-body case. 

Given our understanding of the Hamiltonian structure, canonical quantization is in principle 
straightforward, but the construction of the canonical Hilbert space is not. We have seen that 
in the 2-body case, it is possible to find monodromic cigcnfunctions by projecting over the 
U{1) gauge variable 9 from —00 to +cx), after a weighted sum over monodromies and careful 
regularization. It is also possible that a similar sum and projection (for instance a harmonic 
projection on U{2) or on a proper subgroup) is able to define monodromic wave functions for 
the 3-body system, together with their relevant quantum numbers. Further analysis is needed 
in this direction. 
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A SU{2) monodromies and relabeling symmetry 
A.l Monodromies 

We compute the monodromies of the spinors ctq and in Eq. (j43p and we check that they indeed 
transform according to the matrices M31 and M32 in Eqs. (j45|) and (|46|) when the particles foop 
around each other. While the transformation C, — )• e^™C is straightforward to perform on the 
expressions of f^i, and f2b given in Eq. ([28l) and leads to the diagonal matrix M31 (see Eq. (jlSll ). 
the monodromy around the particle at position ^ = 1 is a bit more tricky. The details being 
quite lengthy, we provide here only the main steps in order for the reader to be able to reproduce 
the full calculation. 

We start from the expressions of f^b and f2b in Eq. ()28p . and we apply the well-known 
hypergeometric transformation [1^ 

^/ , N T(c)T(c — a — b) ^, , 

F{a, b, c, z) = ) -(Fa, b,a + b - c + 1,1 - z) 

i (c — aji [c — 0) 

+ (1 _ zY—bF{c - a, c - 6, 1 + c - a - 6, 1 - z) (133) 

r(a)r(6) 

in order to change the argument of the latter from C, to 1 — C,. Then, since the obtained hyper- 
geometric functions are analytic around the point 1 — C = 0, the monodromy transformations 
may be read from the prefactors. We find 

/Sb = C^^i^^^ -— F /ii, 1 - ^2, 1 - At2 - /^3, 1 - C 

r(^) 

+ ^^^+^3(1 _ ^)M.+^3 r(l + ^3)r(-M2-/.3) ^^^ + ^3, 1 + ^2 + ^3, 1 _ (134) 

As for /2b, a similar transformation may be applied, which leads to 

f r(m + At3)r(i - Ml - /X3)r(^2 + /is) 11 t A^ 

726 = FTTFT^i ^ ^ -/is, 1 - /i, 1 - /i2 - /is, 1 - C) 

r(/^)r(i -/ii) 

, r(/xi + M3)r(M2)r(i - /^i - /^3)r(-//2 - /^s) ^^;,2+^3 ^^..^ 1 , _l t A^ 

+ r(-M3)r(i - /i)r(M) " " ^ + + /^3, 1 - 0- 

(135) 

Then, one applies a further transformation to the two hypergeometric functions which appear 
in the previous expression, 

F(-/i3, 1 - /X, 1 - /X2 - /is, 1 - C) = C'^l+'^^F(1 - ^2, /XI, 1 - /X2 - /is, 1 - C) .^gg. 
F(l - Ml, /i2, 1 + //2 + /is, 1 - C) = C'''+'''i^(/i, 1 + /is, 1 + /i2 + /is, 1 - C) 

which enables one to write /2b in a form that is similar to /sf, as far as the hypergeometric 
functions are concerned: 

_ si r(M2 + /xs)r(//i) 1 „ 1 „ „ 1 n 

726 — (./^l, 1 - /i2, t - /i2 - /is, t - (,j 

Sis r(M) 

_ .12S.Sr(l+Ms)r(-/.2-/is) ^,,^,3(i _ ^)M.+.3^(i + ^3, 1 + /i2 + /.S, 1 - C). (137) 

S1SS2 r(i-/x2) 

One now recognizes that the hypergeometric functions which appear in the transformed expres- 
sions for the /3b and f2b are the same. 
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In this form, the hyper geometric functions are invariant by the monodromy transformation 
— 1 — )■ e'^™{C — !)• Only the prefactors transform, and in a trivial way. Lengthy but straight- 
forward calculations lead to an expression of the relationship between the transformed of the 
vector (/36,/2fe) and its untransformed form through the multiplication by a matrix. Setting 
the relative normalizations in front of f^h and f2b as in Eq. ()43p . we find that this matrix is 
precisely M32 given in Eq. (I45p . Thus ab transforms as a SU{2) spinor (up to a U{1) phase) 
under monodromy transformations. 

A. 2 Label exchange symmetry 

We start from Eq. (144p for the momentum P = {P3,P2), namely 

P = Dz^-^ [at + {Ca - CW] ■ (138) 

The dilation factor D is obviously invariant under any relabeling, see its definition in Eq. (j7]). 

Let us first exchange the labels of particles #1 and #3. This amounts to replacing by 
pi = —p2 — P3, to performing the conformal transformation — )• C/(C — 1) and the substitution 
Z2 — )• 22(1 — 0- The transformed momenta read 

P = Dz^-\l - Cr-^ (^at, + Y^^a) , (139) 

where we have put a "tilde" sign above the transformed quantities. 

We need to understand how cjf, and aa transform. The normalization factors \/T^ and ^/N^ 
are invariant in this case, as seen from their definitions (I32p . As for the transformations of 
the various functions / defined in Eq. (j28p . the key formula is the following identity between 
hypergeometric functions: 

F(a, b, c, C/(C - 1)) = (1 - O'Fic - a, b, c, C). (140) 

Using this formula, expressing aa is then straightforward. Defining 

/gi^Mia o\ ,,,,, 

r3i-( J, (141) 

we obtain 

aa = {I - Cf^'^m ■ CJa. (142) 

As for the transformation of the components /3b and f2b of ai,, one needs to use of the respective 
contiguity relations for hypergeometric functions 

(c - l)F{a, b,c-l,z)- bF{a, b + l,c, z) + {b - c + l)F{a, b, c,z) =0 (143) 

with a = 1 — fi2, b = Hi, c = 1 + fiis, z = and 

cF{a, b, c, z) - bzF{a, b + l,c + l,z) - cF{a - 1, 6, c,z) =0 (144) 

with a = 1 — /i3, 6 = 1 — /u, c=l — /U13, z = C- We arrive at 

ab = {l-0^~''m-i<^b-Caa). (145) 

Combining Eqs. (|142p and (jl45p with Eq. (|139p . one easily finds 

P = T31 • P. (146) 
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We recognize that t^i is the transformation matrix (09]). (Note that the phases stem from the 
replacement of C — 1 by e~*'^(l — C)- The convention with the opposite sign for the phase would 
have led to a different transformation matrix). 

Next, we exchange the labels of particles #2 and ^^3. This amounts to exchanging p2 and 
P3, to transforming into 1/C and to substituting Z2 by Z2C- The transformed momenta thus 
read 



P = Dzi^-\^'-^ l^dk + ^^^^aj • (147) 

The transformations of yJN2 and ^JN^ are straightforward from their definitions (|32|) . Again, 
expressing the various components / of o" with the help of the components of a requires involved 
manipulations of the hypergeometric functions. We first need to use the formula 

F{a, b, c, 1/z) = {e-'^zY F{a, a - c + l,a - b + 1, z) 

r(o)r(c — a) 

+ ^[""f^^^'l (e-^-z)' F(6, b-c+l,b-a + l,z) (148) 
1 (ajl [c — 0) 

in order to express the /'s, which are hypergeometric functions of argument 1/C, as linear 
combinations of the /'s. Without any further transformation, we arrive at 

= C'-'^Tss • CJa, (149) 

with 

The functions f^h, f2b are deduced from the functions /sa, /2a by simply replacing /i2 by /i2 + 1. 
The hypergeometric functions F{—fi2, fj.i, fJ-i3,C) and F(l — ^, 1 — 113,2 — Hi3,C) then appear, 
which are not part of our original basis /a , /2 • But they can actually be expressed as linear 
combinations of the /-functions using contiguity relations. Thanks to the identity 

cF{a, b, c, z) + {a- c)zF{a, b + l,c + 1, z) + (z - l)cF{a, 6 + 1, c) = (151) 

with a = fj.1, b = —fJ-2, c = fiis, z = (, we write 

Fifll, -M2, W3,^) = ^^^'^''^ , C'-^'' [f3b + (1 - C)/3a] . (152) 

r(^i)r(^3) 

Using 

(a - b){a - c + l)F{a, b + l,c,z) + a{a - b){z - l)F{a + 1, 6 + 1, c, z) 

+ ic-l){a-b)F{a,b,c-l,z) = (153) 

with a = 1 — fi, b = —fJ-3, c = 2 — ^13, z = C, we get 

F(l - ^, 1 - ^3, 2 - /ii3, C) = - wT— TT [^26 - (1 - C)/2a] • (154) 

r(l + ^2)r(/il3 - 1) 

Inserting these identities in cjfe expressed with the help of the hypergeometric functions of argu- 
ment we find 

= C'~^T32 ■ K + (1 - Oaa] . (155) 

The replacement of aa and o"b expressed with the help of and ai, (Eq. ()149p and ()155p ) into 
Eq. (I147p leads to the transformation 

P = r32 • P, (156) 

which is an other way to write Eq. ([50]) . 
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